Exact self-accelerating cosmologies in the ghost-free bigravity and massive gravity 
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Within tiie tlieory of tiie gliost-free bigravity, we present tlie most general cosmological solution for 
which the physical metric is homogeneous and isotropic, while the second metric is inhomogeneous. 
The solution includes a matter source and exists for any values of the theory parameters. The 
physical metric describes a universe with the late time acceleration due to the effective cosmological 
term mimicked by the graviton mass. When perturbed, this universe should rest approximately 
homogeneous and isotropic in space regions small compared to the graviton Compton length. In 
the limit where the massless graviton decouples, the solution fulfills the equations of the ghost-free 
massive gravity. 
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Considering theories with massive gravitons [ij is mo- 
tivated by the observation of the current acceleration of 
our universe 0, since the graviton mass can induce an 
effective cosmological term. Although such theories can 
exhibit unphysical features, as for example the Boulware- 
Deser ghost H, the recent discovery of the special mas- 
sive gravity 0| and its bigravity generahzation Q which 
are ghost- free @, suggests that such theories can in- 
deed be good candidates for interpreting the observa- 
tional data. This motivates studying cosmological so- 
lutions with massive gravitons. 

The first self-accelerating cosmologies in the ghost-free 
massive gravity were obtained without matter and de- 
scribe the pure de Sitter universe Q- The matter source 
was then included, but only for special values of the the- 
ory parameters Q . For these solutions the physical met- 
rics is of the Freedman-Robertson- Walker (FRW) type, 
but the fiducial metric is inhomogeneous. This means 
that perturbations around the background can give rise 
to inhomogeneity, although this effect should be sup- 
pressed by the smallness of the graviton mass Q. One 
more similar solutions was found in Ref.(9|, where it was 
argued that solutions of this type should be generic for 
massive gravity. The theory also admits solutions for 
which both metrics are FRW, but these show a nonlinear 
instability and seem to be unphysical LLC 

The generalizations of results of 
ghost-free bigravity were obtained in [11 1 
solutions both metrics are dynamical but are not simul- 
taneously diagonal, and the second metric does not share 
the translational symmetries of the first one. In the bi- 
gravity too, only special solutions of this type are known 
- either witlrout matter or for constrained values of the 
theory parameters. 

In what follows, we present the most general cosmo- 
logical solution for which the physical metric is FRW but 
the second metric is inhomogeneous. We construct this 
solution within the ghost- free bigravity of but it de- 
scribes the massive gravity case as well, since its second 
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metric becomes flat in the limit where the massless gravi- 
ton decouples. The solution includes the matter source, 
it exists for all values of the theory parameters, and its 
physical metric describes a FRW universe that can be 
spatially flat, open or closed, and which shows the late- 
time acceleration due to the effective cosmological term 
mimicked by the graviton mass. 

The ghost-free bigravity— The generic bigravity the- 
ory |13| is defined on a spacetime manifold equipped with 
two metrics g^^i, and whose kinetic terms are chosen 
to be of the standard Einstein-Hilbert form. To describe 
the ghost- free bi grav ity Q, it is convenient to use the 
tetrad formalism |14| , in which the inverse of and the 
/^i^ are parameterized as 
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where R and TZ are the Ricci scalars for g^^^ and /^y, re- 
spectively, G and Q are the corresponding gravitational 
couplings. Sra describes the ordinary matter, which is 
supposed to directly interact only with g^j^^. The interac- 
tion between the two metrics is parameterized as 
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-gd X, 



(3) 



where 



C4 



'■q"/}"7^^5 , (4) 

C3, C4 are parameters, and K^^ = <^i? ~ 7^,y with 

7^. = e>^ (5) 
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These expressions define the theory with two gravitons, 
one of which is massless and another one is massive, with 



2 



the mass = a {1 + G/G). One can introduce the angle 
77 such that the parameters cr, Q are expressed as cr = 
TO^ cos^ 77 and Q ~ Gta.n^ij. Instead of €3,04 one often 
uses m the Uterature = — C3/3 and = — C4/I2. 

Varying the action with respect to and uj^ gives the 
field equations 
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Here and CJ^ are the Einstein tensors for g^^ and /^i/, 
respectively, while 
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These equations should be supplemented by the conserva- 
tion condition for the matter energy-momentum tensor, 
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, where V p is the covariant derivative with 



respect to g^^. 

The field equations imply that T^^^ 
metric, hence so is 7^''' 



g^T^ is sym- 



9^'^l^p; which means that 
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where loa^ — i]AB^p- The latter property implies that 
l^al^u = therefore 7^^ = Vg^'^Uu, which corre- 

sponds to the original theory parameterization QIjS- 

li G then the massless graviton decouples and, if 
only f^i, becomes flat in this limit, the theory reduces to 
the massive gravity of 01 • 

Spherical symmetry— Introducing the spherical coor- 
dinates = (i,r, 7?, (^), the most general expression for 
the two tetrads subject to the condition ([10]) is [ll| 
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where Q, R, a, 6, c, u are functions of t, r. It is straight- 
forward to compute 7''^ in ^ and obtain the following 
non-zero components of in 



Tq = a6 + 2au — 3a + + 04(71 — l)^(a — a6 ^ 
-I- csiu - l){au + 2ab - 3a + 2c^), 



i(u + a + 5- 3) + C4m(1 - u)[c^ + (a - 1)(&- 1)] 



-I- C3u[{a + b - 2)u + + ah - 2a - 2b + 3], 
with rl" obtained from via a -(r^ b, also r,'^ • 
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and 
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The interaction Lagrangian Q reduces to 
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straightforward to evaluate the two energy-momentum 
tensors in ([5]). 

Homogeneity and isotropy— Let us assume the met- 
ric g^i, and the distribution of the ordinary matter to 
be homogeneous and isotropic. This can be achieved 
by setting Q = N = a{t) and R = a(t)/k(r) with 
fkif) = {?■, sin(r), sinh(r)} for k — 0,1,-1, respectively. 
We choose the matter to be of the perfect fluid type with 
8^GT('"\^ = diag[p(t), -P(0, -F(i), -Pit)]. 

Since the Einstein tensor for g^i^ is diagonal, so should 
be the energy-momentum tensor on the right in ([B]), 
therefore one should have = 0, which requires that 

= 0. Now, r" in will vanish if either c = 0, or the 
expression between the brackets vanishes. We shall be 
considering below the case where c ^ and the metrics 
are not simultaneously diagonal, since the c = case has 
already been studied in detail [11| . 

If c 7^ 0, then t° in ^ will vanish if 



C3 + C4 



2C3 + C4 - 1 ± V C§ - C3 + C4 + 1 ) . (15) 



Inserting this into the above formulas, we find that the 
energy-momentum tensors in ([SJ become diagonal, with 
constant 00 and rr components, 

T° = T; = {u- 1)(c3U - u - C3 + 3) = a. 
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To" = T; = -^{C3U - C3 + 2) = A. 
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The Bianchi identities for Eqs.([6]) then imply the conser- 
(s) „ 

vation condition Vp ^ ^ , whose the only non-trivial 
component is (we denote ' = dt and ' = 9^) 
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It is worth noting that a similar condition for Ty fol- 
lows identically, due to the diff.-invariance of S'int , so that 
there is no need to impose it separately. 
Now, using the above formulas one finds 

rrO T-vJ C37J - 7/ - C3 + 2 21 /ISA 

To -T^ = [{u - a){u -b) + c]. (18) 

In view of ((T7| this should vanish, so that either the first 
or the second factor on the right should be zero. The 
former case was considered in Q, [ll| (see also [l6j). 
However, the condition c^u — u — C3 + 2 = constraints 
the possible values of the parameters 03,04, so that the 
solutions obtained in this way are not general. We there- 
fore abandon this condition in what follows and require 
instead that 



{u - a){u - 6) + 



0. 



(19) 
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In view of this, one has r° = T| and = , which 
impUes that both energy-momentum tensors are propor- 
tional to the unit tensor, = Xd^i and = A^^;. The 
field equations (HI),© then reduce to 



K5{ + SttCT^'"^/ , 
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(20) 
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' 77 A. As a result, 



with A = 

the two metrics actually decouple one from another, and 
the graviton mass gives rise to a cosmological term sep- 
arately for each metric. However, one has to remember 
that solutions of ([20]), ((2T|) should in addition fulfill the 
consistency condition (fT9|) . 

Solution for g^i,— Eas. (|20p for g^i, comprise a closed 
system, with no additional conditions imposed, so that 
we can solve them. The metric is 



such that /^^ = rjAB^'^O^ . At the same time, can 
be expanded with respect to u)^ from (|lip . The two sets 
of 1-forms may differ from each other by a local Lorentz 
boost, so that 

= 6*° sec a + 6*^ tan a, cj^ = 6'^ sec a + 6*° tan a, (26) 

where a is the boost parameter. Using the explicit ex- 
pressions for oj^ and 9^ and comparing the coefficients 
in front of dt, dr in (j26p gives four conditions, which de- 
termine a and a, 5, c in terms of A,T,U. As a result, 
the consistency condition (|19p assumes the form 



dsl = a?{t){dt^ - dr^ - fl{r)dn^), 



and the equations reduce to 



A- 



(22) 



(23) 



UT' - TU' - u^a? + uav/A+A_/A = (27) 

with A± = A^T + U' ± {A'^T' + U). This equation de- 
termines T(t, r). 

Let us first consider the 7; — > limit, when A = and 
A = 1, in which case exact solutions of (P7)) can be found. 
For fc = 0. when U — uar, we find 



T{t,r)=C 



*a2 



■dt 



u 
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where p is determined by the matter conservation con- 
dition, p + 3(a/a)(p + P) =0. This describes a uni- 
verse filled with the ordinary matter and containing the 
cosmological term mimicked by the graviton mass. At 
early times the matter density p dominates, but at late 
times the cosmological term wins, which leads to the self- 
acceleration. 

Solution for Eqs. (PT|l should determine the metric 



where C is an integration constant. This solution agrees 
with the one obtained in Q for C3 = C4 = 0, u = 3/2. 
For fc = 1 one has U = ua. sin r and we find 



T{t,r)^ / ^y{C^ +u^){a^ +a^)dt + Cacos{r). (29) 



For fc = — 1 and U = ua sinh r we obtain 



ds} 



[adt + cdrf~a^bdr~cdtf-U^dn\{2A) T{t,r) = ^ {C^ - u^){a? - a^) dt + Cacosh(r). (30) 



Here a, 6, c are free functions of t, r, but U ~ uR{t, r) is 
already fixed by the above considerations. In addition, 
a,b,c should satisfy the constraint (fT9|). One could there- 
fore wonder if the system is not overdetermined and the 
freedom is enough to fulfill all conditions. 

To see that the latter is indeed the case, we notice that 
the function U can be considered as the new radial coor- 
dinate. The temporal coordinate can also be changed, in 
such a way that in new coordinates T, U the metric be- 
comes diagonal. The source term in (|2ip is invariant un- 
der reparameterizations. Therefore, the problem reduces 
to solving the Einstein equations with the cosmological 
constant A to find a diagonal metric parameterized by the 
Schwarzschild coordinate U. The solution is the (anti)de 
Sitter metric 



If A 7^ 1 then exact solutions are difficult to find, how- 
ever, at least when A is small, the solution can be con- 
structed perturbatively as 



(31) 



T = To + ;^(-A/3)"r„. 



Here To corresponds to the zero order expressions ([28])- 
([50)1 , while the corrections T„ can be obtained by separat- 
ing the variables with the ansatz r„ = X^m^o fm(t)g"^{r) 
where g{r) = {r^, cos(r), cosh(r)} for fc = 0,1,-1, 
spectively. For example, for fc = one finds 



re- 
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(32) 



df = A^ dT^ 



dip 



(25) 



where A^ = 1 — AU'^/3. There remains to establish the 
correspondence between the T, U and t, r coordinates and 
to fulfill the constraint ([TO]) . 
Let us introduce 1-forms 

dU 



AdT, 9^ 



9^ 



with C being the same as in ([28]) : similarly for n > 1. 

This completes out constructions, since all field equa- 
tions and the consistency condition are now fulfilled. 
Discussion— We have obtained the cosmological solu- 
tion in the ghost-free bigravity with matter, for arbitrary 
values of the theory parameters (provided that u in (|15p 
is real). One can choose 03,04 such that A > 0. The 
metric g^j, then describes a FRW universe which can be 
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spatially flat, open, or closed. It is matter-dominated 
at early times, but at late times it enters the accelerated 
phase due to the effective cosmological term mimicked by 
the graviton mass. The metric f^^, is the static (anti)de 
Sitter (PSj) parameterized by T, U. 

The two metrics are not simultaneously diagonal and 
do not have the same Killing symmetries. In particu- 
lar, the translational symmetries of g^^ are not shared 
by /^i/, since the Stueckelberg fields 0° = T{t,r), (f)'^ = 
?7(t, r){sin'i9 cos (p, sin t9 sin (y9, cos (/?} are inhomogeneous 
functions of . One can therefore expect the fluctua- 
tions around the background solutions to be inhomoge- 
neous, but this effect will be sourced by the terms pro- 
portional to m? in Eqs.lini),© and so will be small in the 
regions smaller than 1/m, in agreement with the argu- 
ments of 

Setting the matter density to zero, the solution for g^j^y 
is pure de Sitter, rewriting which in the static coordinates 
reproduces the static bigravity solutions found in (l2| . 

The solution exists for any value of ry. When 77 — > 
then A — > 0, so that f^^ becomes flat, while g^y still 
describes the expanding universe. Therefore, we obtain 
in this limit the solution of the ghost-free massive gravity 



with a flat reference metric with the Stueckelberg 
fields expressed by (P5)) - ((50|) . 

Our solution therefore describes the most general cos- 
mology with the FRW g^ u but with inhomogeneous /^y, 
and it equally applies within the bigravity and massive 
gravity contexts. The latter property is quite special, 
since in general the metric /^j^ does not necessary become 
flat for 77 — >■ 0, while generic massive gravity solutions are 
not guaranteed to extend to the bigravity [T5| . 

All other known cosmological solutions can be obtained 
by setting c = in eas. pT|) - p^ . in which case both 
metrics are FRW. The corresponding bigravity solutions 
do not always show the late time acceleration and do 
not cover the massive gravity case , [l^l ■ The massive 
gravity cosmologies with two FRW metrics do not extend 
to the bigravity and show a non- linear instability [loj . 

It seems therefore that our solution is the most sensible 
physically, since it covers all possible cases and shows the 
late time acceleration expected for massive gravitons. 

Note added— When this text was being completed, 
there appeared the article [l^ on massive gravity (77 = 0) 
cosmologies whose analysis is partly similar to the above 
discussion, although it does not give the explicit solution 
for the Stueckelberg scalars. 
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